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ABSTRACT 

We have developed a generalization of the Zeldovich approximation (ZA) that is exact in 
a wide variety of situations, including plannar, spherical and cilyndrical symmetries. We have 
shown that this generalization, that we call complete Zeldovich approximation (CZA), is exact 
to second order at an arbitrary point within any field. For gaussian fields, the third order error 
have been obtained and shown to be very small. For statistical purposes, the CZA leads to 
results exact to the third order. 

Subject headings: cosmology: theory — cosmology: large-scale structure of the Universe - 
gravitation 

Developing a simple analytical approximation that accounts accurately for the non-linear evolution 
of density fields seems to be a rather interesting task. To reconstruct with some accuracy the initial 
conditions from present velocity field in the line of sight one such approximation is needed. To obtain the 
statistical properties of the present field in terms of those for the initial one; the non-linear corrections to the 
microwave background or to the Gun-Peterson effect, an approximation of this kind is highly convenient. 
An approximation that has been widely used to this purposes is Zeldovich approximation (ZA), where 
the density fluctuation, S, is a unique function of the proper values, A*, of the linearly calculated local 
deformation tensor || (u being the peculiar velocity field) . 

(l+sr^na-A,) (!) 

i=i 

This approximation give generaly rather good results and it is particularly convenient for deriving the 
statistical properties of the present field, since we only need the statistical properties of the Xj in the initial 
field. In this approximation the local deformation tensor is that given by the linear theory. So, although it 
is a good approximation to all orders, it is exact only to first order. On the other hand, in the Lagrangian 
perturbative development (LPD) the deformation tensor is formaly exact (Bouchet et al. 1995), but it is 
not a unique function of the A^. Within this context, the question easily arise as to whether is posible 
to find an approximation depending only on the and substancialy more accurate than ZA. Reisnegger 
& Miralda-Escude (1995) considered an extension of ZA (EZA) that is exact for planar, spherical and 
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cilyndrical symmetry. However, although this approximation gives usually better results than ZA, it is 
not fully consistent. Hui and Bertschinger (1996) developed an approximation (LTA) that is exact for any 
fluid element such that the orientation and axis ratios of the gravitational and velocity equipotentials are 
constant along its trayectory. This include planar, spherical and cilyndrical symmetries. The problem with 
this approximation, which gives excelent results, is its complexity, that makes very difficult to determine 
the explicit dependence on the Aj. The purpose of this letter is to present the explicit dependence of 8 as a 
function of the Aj for the most accurate approximation that is a unique function of these quantities. This 
approximation we call complete Zeldovich approximation (CZA). Like EZA and LTA, CZA is exact for 
planar (Ai ^ 0, A2 = A3 = 0), spherical (Ai = A2 = A3) and cilyndrical (Ai = A2, A3 = 0) symmetries, but, 
unlike those approximations, it is also exact in another wide variety of cases that essentialy correspond to 
the collapse of an initial top-hat elliptical density fluctuation. We have shown that CZA is exact to second 
order and computed the third order error, which is very small. We have also shown that the exact evolution 
at an arbitrary point may be expressed in terms of a simple extension of the CZA containing some aditional 
variables (the CZA being recovered when these variables are set equal to zero). All this is shown in detail in 
an accompanying paper; here we shall only present the CZA and comment upon its meaning and derivation. 

Approximations depending only on the A, are usually called local. We shall retain this convention, but 
it must be noted that the Aj, although defined locally, are non- locally generated. The sum of the Ai, which 
is equal to the linear density perturbation, Sl, can take any value at a given point, regardless of the values 
taken at other points. To obtain the Aj, however, we must integrate the continuity equation, so the Ai at a 
given point depends on the whole field 5l {x) . Another way to point out the non-local caracter of the Ai is 
by noticing that the quantities Ai — ^ are generated by the action of the linearly calculated local tidal field 
which, obviously, depends on the whole field 8l{x). Keeping this in mind, it is not surprising that although 
knowing Sl at a point at some initial time let us to obtain the evolution only to first order, knowing the Ai 
let us to obtain the evolution to second order. We shall now describe the steps we have followed to obtain 
the CZA. We choose the ansatz: 

(l + sy^Hii-nWXi), (2) 

i=l 

where Vj(A) are certain functions of the Ai. Within this ansatz, ZA corresponds to the zeroth order 
approximation for ri (ri = 1). It is interesting to note that, although independently developed this ansatz is 
similar to that chosen by the authors of the EZA, except for the fact that they assumed all j-j to be equal; 
an assumption that is incompatible with exactness to second order. 
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To determine the functions r.; we use the constraints imposed on them by: considerations about the 
symmetry of Xi with respect to permutations of the indexes; the fact that for planar symmetry (^) must be 
exact with ri — 1; the form of the exact solution for spherical collapse and compatibility of the form of ri 
with dynamical equations. These conditions determine uniquely. Let us comment them in more detail. 

Rotational invariance imply that the must reduce to each other through permutations of the indexes. 
So, they all derive from the same function, r(u). 



r t (X) = r(u)\ s={XuX . M) (3) 

Furthermore, rotational symmetry in the plane perpendicular to the i-th proper axis imply that Xj, 
Xk must enter symmetricaly in this expression. So, r(u) must be symmetric with respect to its second and 
third arguments. We now assume that a series expansion of r(u) in powers of the Ui exist. We shall see 
later that this series converges for all relevant u. The symmetry considerations we have just mentioned, 
imply that this series can only contain terms of the form: 

oo 

r(u U U2,U3) = 1 + C lm^( u 2+U 3 ) l {u 2 - U 3 ) 2n U™ 

/,m,n— 

p = / + 2n + m, (4) 

where Of are the coefficients of the p-th order terms. Noting that in the planar case (Ai 7^ 0, 

A2 = A3 = 0) expression (|2|) is exact with r< = 1, it is clear that expansion (^) cannot contain terms with 

to ^ and I = n = 0. So, in our notation we must have: 



Hence , there are 2p — 1 terms of order p. 



For spherical collapse, both the actual density fluctuation, S sp , and its linear value, 5l, may be 
expressed exactly as a parametric function of time (Peebles 1980). From these expressions we have derived 
an expression for 6 as an explicit function of 5l' 



x \ - 3 



l + 5 sp = I 1 - r sp (6 L ) — 

r sp (S L ) = 1 + /i(0)y + h(0)^Sl + fa{0)^Sl + / 4 (0)5.86 x 10- 8 *t + / B (fl)2.55 x 1Q- 3 6 5 L + R sp (S L ) 
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Rs P (Sl) = /a(fl)2.58 x 1(T 3 S| ( ^ - 1 j + £ 

V 1 2.065 / 

\E\ < 2 x 10~ 3 for 6 L < 1.57; £ocO(^), (6) 

where 6* stands for all cosmological parameters. For a flat Friedman model all fi(0) are exactly equal 
to one. For a general Friedman model the dependence on £1 (the density in units of the critical one) is very 
mild. When > 1/20 the following is a rather good approximation: 

/;(0) = fi 2 */ 63 (7) 

Comparing expressions (Q) and (j^) and noting that for spherical symmetry Xi = X 2 = A3 = ^r, we 
obtain the following constraint on r(u): 

This relationship imply that the coefficients of order p, Cf m , must satisfy just one equation. So ,for p 
larger than one the coefficients are underdetermined. However, the value of 5 given by (^|) must satisfy the 
dynamical equations (Peebles 1980): 

+ v _ = £ ; 0(f = G a 2 p b (r) / 1 / 1 

at a a J \x — x\ 

->->_•-, dxi(t,qi) 

v = au; u = x; r i (A)A. i = (9) 

oqi 

where Xi are Eulerian comoving coordinates and Xi(t, qi) are the Eulerian coodinates of a particle with 
Lagrangian coordinates q (the last relationship holds in the local proper system). Note that the continuity 
equation is automaticaly satisfied by CZA. From these equations one may see that terms of order p in r 
imply the existence of certain terms of order p + 1. This recursive scheme, together with expresion (^) 
determine completely all coefficients. Their computation, which is not trivial, is given in detail in the 
accompanying paper. Here we simply give the result and comment it. It is interesting to note that the 
process we have followed to determine the coefficients essentialy amounts to analyticaly continue the ?"j(A) 
known in the planar and spherical case in a manner consistent with equations (JoJ) . The situations described 
exactly by the CZA are those where the local deformation tensor (whose proper values are Vj(A)Aj) are 
everywhere the same (or, at least at all fluid elements affecting each other's evolution). This must be so 
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because it is true for the planar and spherical case and preserved by the continuing procedure. For a flat 
Friedman universe, we have found for rj (in the general case terms of order p should be multiplied by f p (0)): 

r t (X t , Xj, X k ) = 1 + -(A, + A fe ) + —(A, + X k f + — (A,- + X k )X t + —(A, - A fe ) 2 + 0.03371(A, + A fc ) 3 

+ 1.63 x \Q- 2 {X,j + X k ) 2 Xi + 2.75 x l(T 2 (Aj + X k )X'f + 10" 3 (Aj - A fc ) 2 A; + 1.2 x 10" 2 (Aj - A fe ) 2 (Aj + A fe ) 

+1.94 x lO" 2 ^ + A fe ) 4 + 9.4 x W~ 3 (Xj + A fe ) 3 A, + 1.58 x lO" 2 ^ + A fc ) 2 A 2 + 1.3 x 10" 2 (Aj + A fe )A 3 
+ 4.3xlO- 3 (A 3 -A fc ) 4 + 8.4xlO- 3 (A J -A fc ) 2 (A,+A A 2 + 7.2xlO- 4 (A,-A A 2 (A,+A A 0A l + i?(A J ^ (10) 

This messy expression includes explicitly terms up to the fourth order although for more purposes 
using up to the quadratic term is enough. Terms of order larger than four are usualy very small. However, 
in the rare cases when they are of some relevance many orders contribute roughly equaly. So, it is not 
convenient to include higher order terms explicitly. Instead, we approximate all those terms by a single 
term, R, that we shall latter give. Expression (|^) with given by ([l(]) give the exact evolution of 5 in a 
field where the local deformation tensor is everywhere given by Ti{X)Xi at a time when the linear one is Ai. 
It is clear that this is the situation within a top-hat cilyndrical fluctuation. So, we may use this case (with 
Q = 1) to check the correctness of the continuing process. In this case we have Ai = 0, A2 = A3 = 4^. So, 
expressions (0) and (|Io|) lead to: 



Sl\ 2 ( 5l S l 



(1 + £cyl) = (^1 -r cy i(£i)y I ; r cy i = r ^y, y,0 
3 107 

= 1 + 28^ + 3528^ + L135 X + 45 X + Rcyl{5L) (U) 

On the other hand, we have obtained r cy i(5i) through direct accurate numerical integration (the error 
of r cy i < 10~ 5 ) and fitted the coefficients in the expansion. These cofficients agree with the predictions 



(those in (11)) well within the fitting errors (0.1% ,0.4% ,3% ,10% for coefficients from the first to the 



fourth). We have also used the numerical results to fit an approximate expression for i? cy i(^L) and find 

i? cy i = 2.2 x lQ- 3 Sf ( 1 - ^] + E 



2.06, 

\E\ < 5 x 10~ 3 (12) 

We may now obtain an expression for R(X) (see expression ([lC])) demanding that the exact result be 
obtained in the planar and spherical cases and that it reduces to a good approximation to i? cy i in the 
cilyndrical case: 
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Aj, Afc) — 

x (V S p(A* + + A fc ) - F sp (A,) + y sp ^ Aj ^ Afc jj + £ 

y sp (x) = r sp (x) - (l + | + + ^x 3 + 5.86 x l(rV) 

= 2.58 x lCrV (l - 1 + E 2 ; E 2 < 2.3 x 1(T 3 for x < 1.57 (13) 

This expression corresponds to = 1 and may be inmediately generalized for an arbitrary cosmological 
model. The maximum error of expression (|l^) with R(X) given by (|l^) is 6 x 10~ 3 being usualy quite 
smaller. For general values of the Ai the situation described exactly by expressions (|^) and (|l0|) is more 
complex than for the three peculier cases already considered. However, for all practical purposes a simple 
generalization of these cases, namely, a top-hat initial ellipsoidal fluctuation, may be considered to be 
described exactly by those expressions. In fact, it may be shown that the intrinsic error of r^(A) in these 
situations as given by ( |l0| ) (»3x 10~ 3 at the time of collapse) is smaller than the error of this expression 
(due only to i?(A)). 

In the top-hat spherical and cilyndrical cases the deformation tensor for outside matter is diferent from 
that for matter within. The same is true for the top-hat eliptical case. However, in this case, unlike in 
the former the outside matter is relevant to the evolution of the matter within. Expression ( [To|) accounts 
exactly to the second order (and a small error to higher orders) for the small contribution of the outside 
matter to the tidal field within the ellipsoid. 

To check the accuracy of the LTA approximation, Hui and Bertschinger (1996) considered the collapse 
of a top-hat initial fluctuation with axial ratios 1:1.25:1.5 and represented (in their fig. 2 ) the evolution of 
the axis predicted by this approximation and that predicted by an approximation (that they called exact) 
that neglects the effect of outside matter. The ellipsoid generates a linear growing mode for the velocity 
field with asociated values given by: Ai = 0.2576a; A2 = 0.3233a; A3 = 0.4191a. (label "1" corresponding 
to the largest axis). As we have said before, expressions (||) and (|Io| ) may be considered exact in this case, 
giving for the evolution of the axis: 



1-9 U 



Aj + A A 



Ai + Aj + A fc 
1.3 



Xi = Wia(l - n(X)Xi) 
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where a is the expansion factor of the universe, and Wi are the axis ratios. This is represented in 
figure 1 along with the predictions of an approximation where all r*j are set equal to their symmetrized 
value ^ ri ( A )+ r2 ( A )+ r3 ( A ) ^ _ This approximation must be very close to EZA (Reisenegger and Miralda-Escude 
1995). It may be shown that by symmetrizing the effect of the outside matter is neglected. This may be 
checked by noting that this approximation is barely distinguishable from "exact" . The LTA and CZA are 
indistinguishable (within a 0.2% ) up to a ~ 1.2, differing very little up to the collapse, that takes place at 
a = 1.584 for CZA and a = 1.613 for LTA. Note that the evolution of the difference between the values of 
the axes given by these approximations is quantitatively the same as that between the exact solution and 
the LTA for an elliptical cloud in empty space (see their fig. 3), such as would be expected if CZA were 
exact. In fact, it may be shown that the error of the value of a at colapse given by CZA is at most 0.3% . 
The agreement between LTA and CZA for a < 1.2 is so complete that if we had an explicit expression (like 
(|Io|)) for the former approximation, the first order coefficient should be very close (within 3%) to that for 
the CZA and those of second order should not differ much. This imply that most likely LTA is exact to 
second order and that it accounts for the effect of outside matter (both facts are related), hence being more 
accurate than "exact" and ECZ (which can not be exact to second order). 



So far we have considered the situations described exactly by CZA. However we are mostly interested 
in the performance of this approximation at an arbitrary point. There is no obvious reason to expect an 
approximation determined by the above considerations to be the best at a random point. To see that this 
is actualy so we first write S at such point in the form: 



r 



.46 x 10- 2 S 2 L Xi 



.34 x 10" 



,.2 1 



10 

y 



(AiA 2 



14 



5 L + 8.46 x 10 _2 <5? 



E 



x i 



(14) 



where the Xi are certain variables defined by the action of some integral operator on Sl (x) and that cannot 
be reduced to functions of the Xi. This expression, that in full would contain more variables, is formaly 
exact, like the LPD (Bouchet et al 1995). But here for each order we have separated the part depending on 
the Xi, that goes into Tj(A)Ai. We may use ( fl4"| ) and the probability distribution of the A^ at points with a 
fixed value of Sl to obtain the statistical properties of S at these points. By comparing what we find with 
the exact results found by Bernardeau (1994) for these statistical quantities we find: 



5>; = 0; (xi) x = C(x 



Sl 
3 



<^>aH^> 2 a + ^(1-CV 
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C = -6.4(C cxp - 0.0544) 



(15) 



The first result is valid at every point and could be derived in other ways ,for example, by comparing 
( |l4| ) with the LPD. The other results, which are of statistical character, give the mean and mean quadratic 
value of Xi over points with a fixed value of the A^. C oxp is a spectral constant defined in the last reference. 

In the context of expression (Q) the CZA, as we have defined it, may be caracterized by the neglect of 
the Xi. This approximation is the same for all fields. However, we could consider a CZA specific for each 
spectrum (for Gaussian fields) simply by inserting in ( fbij ) in the place of Xi , x|, their mean values given in 
( |l5| ) with the value of C oxp corresponding to that spectrum. This approximation give exactly to third order 
the moments of 5 over points with fixed Xi values. Hence, it gives the one point statistics exactly to third 
order. For smooth field (the case considered here) C oxp lies between 0.053 and 0.061 for most interesting 
spectrums. So the general CZA (which is exact to third order for C exp = 0.0544) imply a very small error to 
third order. We have stimated the error of the CZA by computing to third order (the first non-vanishing) 
the RMS fluctuation of the value of 5 over points with fixed Xi values. We found: 



((*- {5)x?)T 



(A? + Xl + Xl)-{X 1 X 2 + ...) — 

7cr 



30, 



3.1 



1/2 



7tJ Sl 



l<-^(C exp -0M71) 

The fact that at every point the sum of the Xi vanishes imply (see ([l4])) that CZA is exact to second 
order. As we have seen, this is most likely to be also true for the LTA, but it is not true for the EZA. 

Both the ZA and the CZA are unique functions of the Xi so, one might wonder where is what makes the 
latter exact to second order. The answer is that in the CZA we use for the proper values of the deformation 
tensor 7*j(A)Aj which is exact to second order, rather than Aj. It must be noted however that the velocity 
field is not given to second order in an explicit manner. The equation V q u = — Y], r%Xj is exact to second 
order but, to obtain the velocity field, we must integrate it. 
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Figure captions: 

Figure 1: Evolution of the axis lengths for a homogeneous ellipsoid embedded in an expanding Universe. 
Initial axial ratios are 1:1.25:1.50. The solid line corresponds to CZA, which is practically exact; and the 
dashed line corresponds to an approximation that makes all equal to their symmetrized value. 
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